Black hole mass and accretion rate of Active Galactic Nuclei 



Xue-Bing Wu, F.K. Liu, R. Wang 
Department of Astronomy, Peking University, Beijing 100871, China 

M.Z. Kong, J.L. Han 
National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China 

(Received February 28 2006) 

The determination of the central black hole mass is crucial to the understanding of AGN physics. 
In this paper we briefly review some methods that are currently used to estimate the black hole 
mass of AGNs. Particularly we demonstrate the importance of two correlations: one between the 
black hole mass and the stellar velocity dispersion and another one between the broad line region 
size and the optical continuum luminosity. Besides applying these relations in deriving black hole 
masses of various types of AGNs, we also employed the fundamental plane of elliptical galaxies to 
estimate the central velocity dispersions of AGN host galaxies, and then the black hole masses of 
AGNs including BL Lac objects. In addition, we derived another empirical relation between the 
BLR size and Hp emission line luminosity from AGNs with the BLR size measured by reverberation 
mapping studies, and argued that more accurate black hole masses of extremely radio-loud AGNs 
could be obtained with it than using the usual R — L • relation because of the jet contribution to 
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the optical continuum. Finally we pointed out that black hole mass estimation is very much helpful 
to determine the accretion rate and understand the accretion process in AGNs. 
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I. INTRODUCTION 

As we know, mass is the most important physics pa- 
rameter of a black hole. Many other properties of the 
black hole are directly related to its mass. Currently 
astrophysical black holes can be divided into three cat- 
egories: the primordial black holes with the mass of 
about 10 15 g existed only in the early universe, the stellar- 
mass black holes with the mass of about 10 solar masses 
(Mq) existed in some X-ray binaries and the supermas- 
sive black holes with the mass of 10 6 to 1O 9 M existed 
in the center of normal and active galaxies. In addition, 
recently there were also some indications of the possi- 
ble existence of the intermediate-mass black holes (with 
the mass of 10 2 to 10 4 M Q ) in some ultra- luminous X- 
ray sources, but this is still under debate. Black holes 
were also suggested to exist in gamma-ray bursts and 
a magnetized accretion disk around the black hole may 
contribute significantly to the energy output of them (Lee 
& Kim 2002; Kim & Lee 2003). In this paper we mainly 
focus on the supermassive black holes in Active Galactic 
Nuclei (AGNs), which consist of many powerful extra- 
galactic sources including quasars, radio galaxies, Seyfert 
galaxies, BL Lacertae objects, LINERs, etc. Supermas- 
sive black holes (hereafter SMBHs) have been suggested 
to exist in the center of AGNs in order to account for 
the huge power of these energetic extragalactic objects 
(Lynden-Bell 1969; Rees 1984). 

In the last two decades, SMBHs have been found in the 
center of about 40 nearby galaxies (including our Milky 
Way) as well, especially with the Hubble Space Telescope 
and other giant ground-based telescopes (See Kormendy 
& Gebhardt 2001 for a recent review). These discoveries 



were made by three dynamical methods, involving the 
studies on the dynamics of stars, gas and extragalactic 
water masers respectively. However, these direct meth- 
ods for black hole mass measurements can not be applied 
to most of AGNs because (1) AGNs are so bright to out- 
shine the stars, and (2) only in several nearby AGNs we 
can observed the circular motions of their central gas and 
water masers. Previously, some indirect methods, such as 
fitting the observed UV/topical 'Big Blue Bump' and the 
X-ray ion K a line in the spectra of AGNs using the ac- 
cretion disk models, have been adopted to estimate the 
central black hole mass of AGNs. But these estimates 
are highly uncertain because they depends on our knowl- 
edge of the disk inclination, accretion rate and radiation 
emissivity law, and all of them are poorly known. Proba- 
bly the most reliable current method for AGN black hole 
mass estimation is the reverberation mapping. Using this 
technique, the size of the broad line region (BLR) can be 
measured using the time lag between the variabilities of 
continuum and emission line fluxes. The black hole mass 
can be then estimated from the BLR size and the char- 
acteristic velocity (determined by the full width at half- 
maximum (FWHM) of the emission line) using the virial 
dynamics. So far, reverberation mapping studies have 
yielded black hole masses of about 20 Seyfert 1 galaxies 
and 17 nearby bright quasars (Wandel et al. 1999; Ho 
1999; Kaspi et al. 2000; Peterson et al. 2004). How- 
ever, even for the reverberation mapping study, there are 
also some uncertainties of the results mainly due to the 
unknown BLR dynamics and geometry (Krolik 2001; Wu 
k Han 2001a). 

With the estimated black hole masses, a tight correla- 
tion between black hole mass and bulge velocity disper- 
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sion (a) has been found for nearby galaxies (Tremaine 
et al. 2002) and for a few Seyfert galaxies as well (Fer- 
rarese et al. 2001). This implies that the Mbh-& relation 
is probably universal for both active and inactive galax- 
ies. Such a tight relation suggests an interesting possibil- 
ity to estimate the central black hole masses for galaxies 
using the measured values of bulge velocity dispersions. 
This straightforward method is particularly important 
for AGNs. Especially BL Lacertae objects, the reverber- 
ation mapping technique cannot be applied because they 
have no or only very weak emission lines in their opti- 
cal spectra. Using the Mbh~& relation may be the only 
possible way to derive the black hole masses for them. 

An empirical relation between the BLR size (R) and 
the optical continuum luminosity at 5100A (^ 5100 jj) nas 
been derived by Kaspi et al. (2000) using the observed 
data of 34 nearby AGNs in the reverberation mapping 
studies. This relation has been frequently adopted to 
estimate the BLR size and then derive the black hole 
masses for large samples of AGNs, including radio-loud 
quasars (Gu, Cao & Jiang 2001; Oshlack, Webster & 
Whiting 2001). However, the optical luminosity of some 
radio-loud AGNs may not be a good indicator of ionizing 
luminosity, which is usually related to the UV/optical ra- 
diation from the accretion disk around the central black 
hole. The powerful jets of blazar-type AGNs may sig- 
nificantly contribute to the optical luminosity. This has 
been confirmed by the discovery of optical jets in AGNs 
by HST (Scarpa et al. 1999), and the detected optical 
synchrotron radiation in many other radio-loud AGNs 
(Whiting, Webster & Francis 2001). Therefore, using 
the R — L ° empirical relation, which was obtained 
based on the sample of mostly radio-quiet AGNs, one 
would significantly overestimate the actual BLR size and 
hence the black hole mass of radio-loud AGNs. Oshlack 
et al. (2002) have shown that their estimated black hole 
masses would be lower if the synchrotron contribution to 
the optical flux is subtracted. 

In this paper, we will present some results in using the 
empirical relations ( Mr w-cr relation and R — L ° re- 
lation) to estimate the SMBH masses of different types of 
AGNs. We will also report some other progress, includ- 
ing using the fundamental plane relation of the elliptical 
host galaxies and a new R — Lh^ empirical relation to 
estimate the black hole masses of radio-loud AGNs. In 
addition, based on the obtained black hole mass we will 
study the accretion rate of some AGNs and discuss the 
accretion physics around black holes. 



II. APPLICATIONS OF THE Mbh-o" RELATION 

A. Black hole masses of Seyfert galaxies and their 
relation with bulge masses 

Nelson & Whittle (1995) have compiled the measure- 
ments of nuclear velocity dispersions of about 70 Seyfert 




5 6 7 8 9 10 

log M BH (M s ) 



FIG. 1: Correlation of V-band absolute bulge magnitudes 
and SMBH masses for quasars, normal galaxies and Seyfert 
galaxies. The solid line represents the least-squares fits to all 
objects. Figure taken from Wu & Han 2001b. 

galaxies. After excluding some LINERs and normal 
galaxies, we got a sample of 33 Seyfert Is and 32 Seyfert 
2s (Wu & Han 2001b). We estimated the SMBH masses 
of these Seyfert galaxies directly using the Mbh-c rela- 
tion given by Gebhardt et al. (2000). The V-band bulge 
absolute magnitude was estimated from the B-band bulge 
absolute magnitude in Nelson & Whittle (1995) by as- 
suming the B-V color of the bulge being the same as the 
total B-V color of Seyfert galaxy. 

The relation between V-band bulge luminosities and 
SMBH masses for 33 Seyfert Is and 32 Seyfert 2s is 
shown in Figure 1. The black hole masses of Seyfert 
galaxies are mostly in the range of 10 6 to 10 8 Af©. In 
Figure 1 we can clearly observe a significant universal re- 
lation between the V-band bulge luminosities and SMBH 
masses for both Seyfert galaxies and quasars, as well as 
for nearby galaxies. For comparison, the SMBH masses 
of 15 quasars and the absolute V magnitudes of their in- 
ner host were taken from Laor (2001). The SMBH masses 
determined by stellar dynamics for 16 nearby galaxies 
were taken from GOO and their absolute bulge V magni- 
tudes were taken from Laor (2001). Adopting the 'typ- 
ical' uncertainties of Af" 9e and Mbh, the least-squares 
fit for all objects gives, 

M bui ge = _ 12 .i4 ± 0.48 - (1.08 ±O.O6)log(M BH /M ), 

(1) 

The Spearman rank-order correlation coefficient is r$ = 
-0.74, which has a probability of P r = 1.3 x 10~ 17 
to occur by chance. Adopting the standard relation 
My i9 ' = 4.83 — 2.51ogLfc„2 ge /L , and the relation be- 
tween the bulge mass and luminosity (Magorrian et al. 
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1998) , from eq. (2) we get, M BH oc M^± 013 . Laor 
(2001) found M BH oc M^±°- 15 for a sample of objects 
including 9 Seyfert galaxies and 15 quasars. Our re- 
sult shows that this nonlinear relation is more evident 
from a larger sample including more Seyfert galaxies. 
Using another well-selected sample including 22 Seyfert 
Is and 15 Seyfert 2s (Wu & Han 2001b), we also found 
Mbh cx A^jge ±0 ' 14 ^- Therefore, the nonlinear Mbh- 
Mbuige relation is confirmed by our studies. 

B. Fundamental plane relation and black hole 
masses of BL Lac objects, radio galaxies and quasars 

AGNs usually have very bright nuclear emission, which 
makes it difficult to measure their stellar velocity disper- 
sions with the spectroscopic method. It is well known 
for ellipticals that three observables, namely the effective 
radius, the corresponding average surface brightness and 
the central velocity dispersion, follow a surprisingly tight 
relation (so-called the fundamental plane). Some stud- 
ies have shown that the elliptical hosts of radio galaxies 
follow the same fundamental plane as normal ellipticals 
(Bettoni et al. 2001). For about 300 normal ellipticals 
and radio galaxies, Bettoni et al. (2001) found that the 
fundamental plane can be robustly described as 

logi? e = (1.27 ±0.04) log a + (0.326 ±0.007) < fi e > R 

-8.56 ±0.06, (2) 

where R e is the effective radius in kpc, a is the central ve- 
locity dispersion in fcms -1 , and < u, e >r is the average 
surface brightness in R-band. Because the fundamen- 
tal plane probably exists also for the elliptical hosts of 
AGNs, this provides a possible way to estimate the cen- 
tral velocity dispersions and then the SMBH masses of 
AGNs. 

The BL Lac snapshot survey using the HST has ob- 
tained images for 110 BL Lac objects (Scarpa et al. 

1999) . Urry et al. (2000) has obtained host galaxy pa- 
rameters for 72 of these BL lacs. Among them, 63 ob- 
jects have measured redshifts. 51 of them are classified 
as high-frequency peaked BL Lacs (HBL) and 12 of them 
as low- frequency peaked BL Lacs (LBL). Figure 2 shows 
the histograms of our derived black hole masses based on 
the fundamental plane and the Mgff-ij relation for HBLs 
and LBLs (Wu, Liu & Zhang 2002). It is clear that there 
is no significant difference in the black hole masses be- 
tween HBLs and LBLs. Most BL Lacs have black hole 
masses in the range of 1O 7 5 M0 to 10 9 ' 5 M Q . 

A deep HST imaging study of the host galaxies of a 
sample of 10 radio galaxies (RGs), 10 radio-loud quasars 
(RLQs) and 13 radio-quiet quasars (RQQs) has been per- 
formed (McLure et al. 1999). The hosts of both radio- 
loud AGNs and bright radio-quiet AGNs were found to 
be virtually all massive ellipticals. Therefore, we can 
also derive the central velocity dispersions and black hole 
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FIG. 2: Histogram of the derived SMBH mass distribution of 
HBLs and LBLs using the fundamental plane and the Mbh-ct 
relation. Figure taken from Wu et al. (2002). 
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FIG. 3: Histogram of the derived SMBH mass distribution of 
radio galaxies, radio-loud and radio-quiet quasars using the 
fundamental plane and the Mbh-ct relation. Figure taken 
from Wu et al. (2002). 



masses for these AGNs based on the morphology param- 
eters of their host galaxies (Wu et al. 2002). Figure 3 
shows the histograms of the derived black hole masses 
for these AGNs. It is clear that there are no significant 
differences in the SMBH masses among RGs, RLQs and 
RQQs. Most of these AGNs have black hole masses in 
the range of 10 7 5 M Q to 1O 9 M . Our results indicate 
that there is no difference in SMBH masses of BL Lacs, 
RGs and RLQs. Different from some previous claims 
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that RLQs have more massive SMBHs than RQQs (Laor 
2000), our result shows that there is only a weak differ- 
ence in our derived SMBH masses for RQQs and RLQs 
in this sample if their host galaxies are ellipticals. 

Using the same method, we also estimate the primary 
black hole mass of another well-known BL Lac object 
OJ 287 (Liu & Wu 2002). Our result gives the mass 
about 4 x 10 8 M Q , which is consistent with the upper 
limit (10 9 M o ) obtained by Valtaoja et al. (2000) based 
on a new binary black hole model for OJ 287. 



III. APPLICATION OF THE R 
RELATION 
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Reverberation mapping studies of 34 AGNs revealed 
a significant correlation between the broad line region 
(BLR) size and the optical continuum luminosity (Kaspi 
et al. 2000), namely, 
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This empirical relation provides a possible way to esti- 
mate the black hole mass of an AGN from a single optical 
spectrum, from which the continuum luminosity can be 
easily measured. On the other hand, the width (FWHM) 
of Balmer lines, which can be used to estimate the veloc- 
ity dispersion of the BLR. Therefore, The virial mass of 
the central black hole can be estimated with 



M BH = 1-464 x 10 5 ( 
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where Vfwhm is the FWHM value of Balmer emission 
lines. 

Here we show some of our results in estimating the 
black hole mass of double-peaked broad line AGNs using 
the R — L = relation (Wu & Liu 2004). These AGNs 

5100.4 v ' 

represent only a small fraction of AGNs and the double- 
peaked emission lines are usually attributed to the line 
emission from accretion disks. So far about 150 double- 
peaked broad emission line AGNs have been discovered. 
In a completed survey, Eracleous & Halpern (1994, 2003) 
have found 26 radio-loud double-peaked broad emission 
line AGNs. Recently Strateva et al. (2003) presented 
a new sample of 116 double-peaked Balmer line AGNs 
selected from the Sloan Digital Sky Survey (SDSS). In 
addition, double-peaked broad emission lines have been 
also found in some nearby galaxies, including NGC 1097, 
M81, NGC 4450, and NGC 4203. These objects are all 
classified spectroscopically as Type 1 LINERs with very 
smaller nuclear luminosity (< 10 43 erg/s) and very lower 
Eddington ratios (< 10~ 3 )(Ho et al. 2000). 

Strateva et al. (2003) have given the magnitudes 
(Galactic extinction corrected), redshifts and FWHM 
values of the double-peaked H a line for 109 SDSS AGNs 
(7 objects have no FWHM values among 116 AGNs 
in their table 3). From the g and r magnitudes of 
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FIG. 4: The histograms of the black hole masses and Edding- 
ton ratios of double-peaked AGNs in the SDSS and radio-loud 
AGN samples. Figure taken from Wu & Liu (2004). 



these AGNs (available in Tables 1 and 2 in Strateva 
et al. (2003)), we can estimate the rest frame lumi- 
nosity at 5100A Using the conversion formula between 
the FWHM values of H Q and H a \mes,VFWHAi(Ha) = 
.873Vfw h m (Hp) , we can estimate the black hole masses 
of the double-peaked SDSS AGNs with Eqs (3) and (4). 
The same method can be applied to 26 double-peaked 
AGNs found from a survey of radio-loud emission line 
AGNs (Eracleous & Halpern 1994, 2003). The Eddington 
ratios [Lboi/LEdd, where L^oi is the bolometric luminos- 
ity and LEdd is the Eddington luminosity) of these AGNs 
can be also derived if we adopted Lf, i ~ ^^ 51Q0 ^ (Kaspi 
et al. 2000). For 109 SDSS double-peaked AGNs, the 
black hole masses of are in the range from 4 x 10 7 M Q to 
5 x 10 9 M Q , and the Eddington ratios are between 0.002 
and 0.2. The average black hole mass of these AGNs is 
10 8 ' 74 Mq and the average Eddington ratio is 10~ 182 . For 
the 26 radio-loud double-peaked AGNs, the black hole 
masses are from 3 x 1O 7 M to 3 x 10 9 M Q and the Edding- 
ton ratios are between 0.001 and 0.08, The average black 
hole mass of these 26 double-peaked AGNs is 10 8 78 M Q 
and the average Eddington ratio is 10 -2 01 , both being 
similar as those obtained for the double-peaked SDSS 
AGNs. Figure 4 summarized our results. Clearly, our 
derived Eddington ratios for the SDSS and radio-loud 
double-peaked AGNs are larger than the LINER-type 
double-peaked sources, which may reflect the difference 
in the accretion process for these AGNs (see section V 
for discussions). 
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FIG. 5: The upper panel shows the R — L Hfj relation. The 
line shows the OLS bisector fit to the data. The lower panel 
shows the R — L ° relation, which is identical to that in 

5100A 

Kaspi et al. (2000) except that we include Mrk 279 (open 
triangle). The line represents the linear fit given in Kaspi et 
al. (2000). The open and filled symbols denote Seyferts and 
quasars respectively. Figure taken from Wu et al. (2004). 



IV. BLACK HOLE MASSES DERIVED FROM 
THE BLR SIZE AND H,3 LUMINOSITY 
RELATION 

A. The R — Lh relation 

Using the available data of BLR sizes and Hp fluxes for 
34 AGNs in the reverberation mapping studies, we can 
investigate the relation between the BLR size and the 
emission line luminosity (Wu et al. 2004). The Hp flux 
data were available for 16 PG quasars and 17 Seyfert 1 
galaxies. Because there is no available data of Hp flux for 
PG 1351+640 (Kaspi et al. 2000), we exclude this ob- 
ject from our investigation. In addition, we add another 
Seyfert galaxy Mrk 279 in our sample because both the 
BLR size and the flux have been measured recently 
(Santos-Lleo et al. 2001). All Hp luminosity data have 
been corrected for Galactic extinction using the values 
from NED. With these data, we derive an empirical re- 
lation between the BLR size and luminosity: 

Log R (It - days) = (1.381 ± 0.080) + 

(0.684 ±0.106)£o 5 (L^/10 42 ergs s" 1 ). (5) 

The Spearman's rank correlation coefficient of our R — 
Lh relation is 0.91, slightly higher than 0.83 for the 
R — L - relation (Kaspi et al. 2000). 
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FIG. 6: Upper panel: Comparison of the black hole masses of 
radio-loud quasars estimated with R — Lh„ and R — L 

H f 5100A 

relations. The squares represent 59 quasars in Brotherton 
(1996) and the triangles represent 27 PHFS quasars in Osh- 
lack et al. (2002). The diagonal line shows the results are 
identical. Lower panel: The ratios of black hole masses es- 
timated with two different relations are plotted against the 
radio-loudness of radio-loud quasars. The dashed line shows 
that two black hole masses are identical. Figure taken from 
Wu et al. (2004). 



In Fig. 5 we show the dependence of the BLR size 
on L u a and L ° . Obviously these two relations are 

_ 5100A J 

similar, which means that the R — Lh^ relation can be 
an alternative of the R — L ° relation in estimating 
the BLR size for radio-quiet AGNs. 



B. Comparison of black hole mass estimation of 
AGNs from two R-L relations 

Since the R — L ° relation obtained by Kaspi et al. 

5100 A J l 

(2000) has been frequently used to estimate the BLR size 
and the black hole mass of both radio-quiet and radio- 
loud AGNs, it is important to investigate the applicabil- 
ity of such an approach for radio-loud objects. 

Brotherton (1996) studied the emission line properties 
of 59 radio-loud quasars. We adopted his published val- 
ues of absolute V-band magnitude, equivalent width and 
FWHM of Hp emission line. The continuum luminos- 
ity at 5 100 A and the Hp luminosity were calculated af- 
ter considering the Galactic extinction and K-correction 
(optical spectral index was assumed to be 0.3). We 
then estimated the BLR size using both R — Lh and 
R— L ° relations and derived the black hole mass with 

5100A 
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the formula M BH = 3V£ WHM R/4G (h ere we assume the 
BLR velocity V ~ (y/3/2) Vfwhm as did by Kaspi ct 
al. 2000). With these two relations, we also estimated 
the black hole masses of another 27 radio-loud quasars 
with available data of both the equivalent width and 
FWHM of Hp emission line in the Parkes Half-Jansky 
flat-spectrum Sample (PHFS) (Drinkwater et al. 1997; 
Francis, Whiting & Webster 2000; Oshlack et al. 2001). 
We compared the black hole masses obtained with the 
R— Lh„ and R— ^ 5100 ^ relations in Figure 6. Evidently 
the masses obtained with the R — Lh p relation are sys- 
tematically lower that those obtained with the R—L 

J 5100A 

relation for some extremely radio-loud quasars. In Figure 
6 we also show how the difference of black hole masses 
varies with the radio-loudness for these two samples of 
radio-loud quasars. It is clear that the difference of black 
hole masses is smaller when the radio-loudness is small 
but becomes larger as the radio-loudness increases. For 
some individual quasars with higher radio-loudness, the 
black hole mass estimated with the R — L ° relation 

5100A 

can be 3~10 times larger than that estimated with the 
R — Lh p relation. For radio-quiet AGNs, however, both 
the optical continuum and emission line luminosities are 
probably free from the jet contributions and therefore 
both can be good tracers of photo-ionization luminos- 
ity. Using these two relations will give almost identical 
results. 

More recently Kong et al. (2006) extended such a 
study to the broad UV emission lines Mgll and CIV. 
Their results confirm again that the relation between the 
BLR size and broad emission line luminosity can be used 
to derive more accurate black hole mass for extremely 
radio-loud AGNs. 



V. DISCUSSIONS 

We have presented some applications of using the em- 
pirical Mbh — cr and R— L relations to estimate the black 
hole mass of various types of AGNs and reported some 
progress in improving these studies by using the funda- 
mental plane of the elliptical host galaxies of AGNs and 
the relation between the BLR size and broad emission 
line luminosity. Although there are still some uncertain- 
ties, these multiple methods can give consistent results 
on the black hole mass of most AGNs and have been used 
to reliably estimate the black hole mass of larger samples 
of AGNs. 

Knowing black hole mass is important to study AGN 
physics. One of such examples is to use the black hole 
mass to estimate the Eddington ratio (Lt> i/ Lsdd, which 
usually defines the dimensionless accretion rate of AGNs. 
From the accretion disk physics, we know that at differ- 
ent accretion rate, the accretion flow may have different 
accretion modes (Abramowicz et al. 1995). Therefore, 
different Eddington ratios of AGNs may reflect different 



accretion physics. Fig. 7 shows the dependence of the 
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FIG. 7: Relation between the broad band spectral index 
aox and Eddington ratios for three different samples (SDSS, 
Radio-loud and LINER-type) of double-peaked broad line 
AGNs (see section III). Figure taken from Wu & Liu (2004). 



broad band spectral index aox on the Eddington ra- 
tios for double peaked broad line AGNs with different 
luminosity. Clearly, the AGNs at higher accretion rate 
have different spectral shape in the overall spectral en- 
ergy distributions (SED) from those at lower accretion 
rate. Wu & Liu (2004) has demonstrated that the ob- 
jects with higher aox values display big blue bumps in 
the optical/UV bands, which is widely believed to be 
from the thermal radiation of an optically thick accre- 
tion disk (Shakura & Sunyaev 1973). On the contrary, 
the objects with smaller aox values (usually with smaller 
Eddington ratios) may host an optically thin, radiatively 
inefficient advection dominated accretion flow (Narayan 
& Yi 1994). Therefore, our knowledge of the accretion 
rate based on the black hole mass estimates is crucial to 
understand the accretion physics. 
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